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Thermoelectric  materials  are  promising  in  converting  heat  directly  into  electricity,  and 
composite  materials  are  often  used  to  enhance  thermoelectric  figure  of  merit.  In  this 
work,  we  develop  a  nonlinear  asymptotic  homogenization  theory  to  analyze  the  effective 
behavior  of  layered  thermoelectric  composite  with  coupled  transport  of  electricity  and 
heat.  The  nonlinearly  coupled  thermoelectric  transport  equations  are  homogenized  using 
asymptotic  analysis,  from  which  the  macroscopic  field  distributions  are  derived  with  local 
fluctuation  averaged  out,  and  overall  thermoelectric  conversion  efficiency  is  established 
using  an  idealized  thermoelectric  module.  It  is  discovered  that  the  thermoelectric  field 
distributions  in  the  composite  are  different  from  those  in  a  homogeneous  material,  and 
they  are  difficult  to  be  fitted  by  homogeneous  solution.  Furthermore,  it  is  noted  that  while 
the  effective  thermoelectric  properties  of  the  composite  can  be  defined  through  a  set  of 
equivalency  principle,  these  effective  properties  depend  on  specific  boundary  conditions, 
resulting  in  effective  figure  of  merit  that  is  not  correlated  with  thermoelectric  conversion 
efficiency  directly.  The  analysis  thus  sheds  considerable  insight  into  the  effective  behavior 
of  thermoelectric  composites  for  their  design  and  optimization. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 


Thermoelectric  materials  are  promising  in  converting  heat  directly  into  electricity  and  vice  versa,  and  have  attracted 
significant  interests  in  recent  years  for  their  potential  applications  in  waste  heat  recovery  (Yang  and  Caillat,  2006;  Narducci, 
2011),  solid  state  cooling  and  thermal  management  (Tritt  and  Subramanian,  2006),  solar  energy  harvesting  (Kraemer  et  al„ 
2011),  and  carbon  reduction  (Bell,  2008).  The  efficiency  of  thermoelectric  conversion  is  governed  by  a  dimensionless  figure 
of  merit  ZT  (Tritt  et  al.,  2008): 


0) 


where  a,  a,  and  k  are  Seebeck  coefficient,  electric  conductivity,  and  thermal  conductivity,  respectively,  and  T  is  the 
temperature.  High  thermoelectric  conversion  efficiency  requires  not  only  high  Seebeck  coefficient,  but  also  high  electric 
conductivity  and  low  thermal  conductivity,  and  this  turns  out  to  be  rather  difficult  to  achieve  simultaneously,  since  all  these 
properties  are  intimately  related  to  each  other.  For  example,  a  good  electric  conductor  is  usually  a  good  thermal  conductor 
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as  well.  As  a  result,  one  of  the  main  strategies  in  developing  high  performance  thermoelectrics  is  to  engineer  hybrid 
composites  to  enhance  their  conversion  efficiency,  especially  with  nanostructure  engineering  to  take  advantages  of  size  and 
interface  effects,  and  large  enhancement  in  thermoelectric  Figure  of  merit  has  indeed  been  reported  for  a  number  of 
composite  systems  (Poudel  et  al„  2008;  Heremans  et  al.,  2002;  Poudeu  et  al.,  2010;  Pei  et  al.,  2011;  Xie  et  al.,  2009,  2010; 
Zhou  et  al.,  2008;  Zhu  et  al„  2009;  Kim  et  al.,  2007;  Lin  et  al.,  2005;  Gothard  et  al„  2008;  Alboni  et  al.,  2008;  Cao  et  al.,  2008; 
Arachchige  et  al.,  2008;  Ke  et  al„  2009;  Vashaee  and  Shakouri,  2004). 

While  vast  amount  of  experimental  works  focus  on  developing  thermoelectric  composites,  there  have  been  only  very 
limited  theoretical  efforts  toward  the  continuum  analysis  of  their  effective  behavior  (Liu,  2012;  Hao  et  al.,  2012;  Yang  et  al., 
2012,  2013).  Using  variational  principle  similar  to  Hashin-Shtrikman  bound  in  elasticity  (Hashin  and  Shtrikman,  1962),  it 
was  argued  that  the  effective  figure  of  merit  of  a  composite  is  bounded  by  the  figure  of  merit  of  its  constituents  (Bergman 
and  Levy,  1991 ;  Bergman  and  Fel,  1999),  though  the  analysis  was  built  on  linearized  thermoelectric  transport  equations,  and 
thus  the  conclusion  needs  to  be  further  examined;  we  will  come  back  to  this  point  later.  Furthermore,  heat  flux  was  often 
assumed  to  be  divergence  free  (Webman  et  al.,  1977),  which  is  also  not  appropriate  for  thermoelectric  materials  with 
coupled  transport  of  electricity  and  heat.  These  motivate  us  to  examine  the  effective  behavior  of  thermoelectric  composites 
using  more  realistic  nonlinear  analysis. 

As  we  discussed,  the  difficulty  in  continuum  analysis  of  thermoelectric  is  the  nonlinearly  coupled  electric  conduction  and 
heat  transfer,  which  have  not  been  sufficiently  elucidated  in  the  literature.  In  fact,  “there  is  a  lot  of  confusion,  particularly  in 
the  composite  materials  community,  as  to  the  appropriate  form  of  the  thermoelectric  equations”  (Milton,  2002).  Recently, 
we  have  carried  out  rigorous  nonlinear  analysis  of  thermoelectric  transport  in  bilayered  composite  structure,  where  we 
showed  that  thermoelectric  conversion  efficiency  higher  than  both  of  the  constituents  can  be  achieved  (Yang  et  al.,  2013),  in 
contrast  to  the  previous  claim.  The  bilayered  composite,  however,  is  a  structure  without  macroscopic  homogeneity,  and  thus 
in  this  work,  we  analyze  the  thermoelectric  composite  with  periodic  microstructure,  as  schematically  shown  in  Fig.  1.  In 
Section  2,  the  general  framework  of  thermoelectricity  is  deduced  from  irreversible  thermodynamics,  and  one-dimensional 
field  analysis  is  presented  for  homogeneous  materials.  Asymptotic  homogenization  (Bakhvalov  and  Panasenko,  1989; 
Kalamkarov  et  al.,  2009;  Triantafyllidis  and  Bardenhagen,  1996;  Homentcovschi  and  Dascalu,  2000;  Tarn,  1997)  is  then 
carried  out  in  Section  3,  and  governing  equations  of  coupled  thermoelectric  transport  are  homogenized.  This  enables  us  to 
analyze  the  effective  behavior  and  conversion  efficiency  of  the  composite  in  Section  4,  with  numerical  results  and 
discussions  presented  in  Section  5.  Key  observations  from  the  analysis  include  the  drastically  different  governing  equations 
and  functional  distribution  of  thermoelectric  fields  in  the  composite  from  those  of  homogeneous  materials,  and  the  ill- 
defined  composite  figure  of  merit  that  is  found  to  be  uncorrelated  with  the  thermoelectric  conversion  efficiency.  Such 
analysis  makes  it  possible  to  optimize  the  conversion  efficiency  of  composite  materials,  which  is  currently  undergoing. 


2.  Thermoelectric  effects 


2.1.  General  framework 


The  continuum  theory  of  thermoelectrics  can  be  traced  back  to  thermodynamics,  where  the  fundamental  relation  defines 
S,  the  entropy  of  a  system,  in  terms  of  its  extensive  parameters  internal  energy  U,  volume  V,  and  the  mole  numbers  of  the 
chemical  components  N,-: 

S=S(U,V,Ni).  (2) 
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We  restrict  ourselves  to  volume  preserving  processes,  and  consider  electron  as  the  only  type  of  charge  carrier,  so  that  the 
fundamental  relation  can  be  recast  as 

s=s(u,n),  (3) 

where  s,  u,  and  n  are  entropy,  energy,  and  mole  number  of  electrons  per  unit  volume,  respectively.  By  differentiating  this 
specific  fundamental  relation,  we  derive 

ds=^du-~dn,  (4) 

where  7  is  the  absolute  temperature,  and  p  is  the  electrochemical  potential  of  electrons.  We  are  particularly  interested  in  the 
flows  of  energy  and  electrons,  for  which  the  rate  of  entropy  production  is  evaluated  as 


Sfell-Jtf-Vj-Jtt,  (5) 

where  Ju  and  Jw  are  the  fluxes  of  energy  and  electrons,  respectively,  while  Vl /T  and  -Vp/T  are  corresponding  affinities  that 
drive  such  flows.  For  the  Markoffian  system  (Callen,  1960)  where  the  fluxes  depend  only  on  the  instantaneous  affinities,  it 
was  proposed  that 

-J^L'nV^  +  L'izV;*, 

J  (/  =  L'12V^  +  L'22vI,  (6) 

where  isotropic  symmetry  is  assumed  for  simplicity,  and  L'y  are  transport  coefficients  with  L'y  =  L  j,-  due  to  the  Onsager 
theorem  (Callen,  1960).  This  set  of  equation  is  equivalent  to  Eq.  (2.17)  in  Milton  (2002),  and  it  is  evident  that  fluxes  of  energy 
and  electrons  are  coupled  by  L'\2- 

For  the  system  we  are  considering,  the  energy  is  transported  through  both  heat  and  electrochemical  potential,  and  thus 

Ju  —  Jq  +  kIiVj  (7) 

where  JQ  is  the  heat  flux.  This  can  be  used  to  rewrite  the  transport  equations  as 
-JN  =  L,iiv^ta2v|, 

Jq  =  f-i2  +  L22V  (8) 

L„=L'„,  Ln  =  pL'u+L'u,  L22=p2L'u+2pL\2+L'22. 

Assume  that  the  chemical  potential  of  electrons  is  independent  of  temperature  such  that 

P  =  ed>,  (9) 

where  e  is  the  charge  of  electron  and  <f>  is  the  electric  potential,  and  notice  that  the  electric  current  density  J  is  related  to 
electron  flux  JN  as 

J  =  ejjy,  (10) 

we  derive  the  coupled  transport  equations  for  electric  current  density  and  heat  flux 


Jq  = 


eli  2 

~T~ 


V^-^2  VT, 


(H) 


and  these  are  equivalent  to  similar  equations  in  Milton  (2002)  as  well.  This  set  of  transport  equations  allows  us  to  establish 
connections  between  transport  coefficients  Ly  with  more  familiar  thermoelectric  properties  (Goupil  et  al„  2011),  including 
electric  conductivity  a,  thermal  conductivity  k,  and  Seebeck  coefficient  a.  In  particular,  notice  that  electric  conductivity  is 
measured  under  isothermal  condition,  such  that 


JUn=o  =  e*Lu 

~ '  V#  ”  T  ' 

Similarly,  thermal  conductivity  is  measured  under  open-circuit  condition,  resulting  in 


(12) 


Jo}j=t0  fltf-23-f-i2 

"  vr'=  * 


(13) 
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Furthermore,  Seebeck  coefficient  is  measured 
condition,  leading  to 


V0l;  =  o_  L 
VT  ~  eTL, 


electric  field  induced  by 


temperature  gradient,  under  open-circuit 
(14) 


With  these  connections,  the  thermoelectric  transport  equations  can  then  be  recast  into  more  familiar  form  as 


-J  =  crVi j>  +  aaVT, 

Jq  =  —  TaaVtp— (ToPo  +  k)VT  =  Taj— xVT,  (15) 


which  are  nonlinearly  coupled  due  to  the  terms  TVrp  and  TVT.  In  the  absence  of  thermoelectric  effect  where  a  =  0,  the 
uncoupled  transport  equations  of  electricity  and  heat  are  recovered. 

It  is  also  evident  from  energy  transport  that 

Ju=Jq  +  0J.  (16) 

and  for  a  system  wherein  both  charges  and  energy  are  conserved,  current  density  and  energy  flux  are  divergence-free, 

V  J  =  0, 

VJ„  =  0,  (17) 


while  the  heat  flux  is  no  longer  divergence-free  in  general, 

V  -Jq  =  -V  ■  (<p\)  =  — •  J,  (18) 

where  Joule  heat  serves  as  heat  source.  This  is  in  contrast  to  a  normal  heat  transfer  problem  uncoupled  from  electric 
conduction. 

Before  we  proceed,  we  make  the  following  observations  on  various  forms  of  constitutive  equations,  (6),  (8),  (11 )  and  (15). 
They  of  course  are  all  equivalent,  yet  subtle  difference  exists.  The  first  set  of  equations,  (6),  appears  to  be  linear,  if  the 
involved  kinetic  coefficients  L'y  can  be  viewed  as  constants  independent  of  temperature  and  electrochemical  potential; 
these  would  make  homogenization  straightforward  (Milgrom  and  Shtrikman,  1989;  Milton,  2002).  Indeed,  as  recently 
pointed  out  by  Liu  (2012),  these  kinetic  coefficients  can  be  regarded  as  constants  under  small  deviation  of  the  temperature 
and  electrochemical  potential  from  their  equilibrium  values,  with  which  the  transport  equations  can  be  linearized.  However, 
under  more  general  circumstances,  these  kinetic  coefficients  are  functions  of  local  intensive  parameters  (Callen,  1960),  and 
thus  cannot  be  viewed  as  constants,  making  the  governing  equation  intrinsically  nonlinear.  These  can  be  better  appreciated 
from  Eq.  (12).  If  L'n  =  Ln  =  oT /e2  would  be  a  constant  independent  of  temperature,  then  the  electric  conductivity  a  would 
vary  inversely  with  respect  to  temperature,  as  noted  by  Liu  (2012).  For  semiconductors  that  most  good  thermoelectrics 
belong  to,  it  is  well  known  that  the  intrinsic  electric  conductivity  generally  increases  with  respect  to  temperature.  In  fact, 
electric  conductivity,  thermal  conductivity,  and  Seebeck  coefficient  all  have  complicated  variation  with  respect  to 
temperature,  and  thus  thermoelectric  transport  equations  have  to  be  regarded  as  nonlinear  in  general.  Eq.  (6)  also  appears 
to  be  advantageous  compared  to  the  other  three  sets  of  equations  in  that  both  JN  and  are  divergence-free,  yet  the  physical 
meanings  of  the  corresponding  kinetic  coefficients  are  less  straightforward.  In  this  work,  we  choose  to  work  with  Eq.  (15), 
for  its  close  connection  with  widely  used  thermoelectric  properties.  It  is  worth  noting  that  Eq.  (15)  is  essentially  identical  to 
the  transport  equations  used  by  Bergman  and  Levy  (1991)  and  Bergman  and  Fel  (1999),  who  linearized  the  equations  by 
treating  k/T  +  era2  as  a  constant,  making  it  possible  to  use  Hashin-Shtrikman  variational  principle  to  bound  the  effective 
figure  of  merit  of  thermoelectric  composite.  As  pointed  out  by  Liu  (2012),  this  would  also  require  small  deviation  of 
temperature  and  electrochemical  potential  from  their  equilibrium  values.  Furthermore,  the  linearized  version  of  Eqs.  (6)  and 
(15)  are  no  longer  equivalent,  and  their  solutions  to  specific  problems  are  different.  In  this  work,  we  do  not  invoke  these 
assumptions  for  linearization,  and  will  examine  the  implication  of  nonlinearity  in  the  transport  equations  to  the 
thermoelectric  properties  of  composites. 


2.2.  One-dimensional  system 


To  be  specific,  we  confine  our  scope  to  one-dimensional  (ID)  system  in  this  work,  wherein  all  the  field  variables  and 
materials  parameters  are  assumed  to  be  dependent  only  on  coordinate  x,  and  independent  ofy  and  z,  and  the  thermoelectric 
transport  equations  are  simplified  as 


?L 

J  a  dx  +  aa  dx  ’ 


(19) 


h  =  «V-*§'  (20) 

In  addition,  the  divergence-free  of  current  density  implies  that  J  =  Ul  =  const  in  ID,  and  the  field  equation  for  heat  flux  is 
simplified  as 

4/q _ tty, 

dx  ~  dxJ' 


(21) 
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The  transport  equations  can  be  combined  with  the  field  equation,  resulting  in  the  following  governing  equations  for 
temperature  and  electric  potential  in  a  homogeneous  medium  wherein  all  the  thermoelectric  properties  are  constant, 


dx1  ™  ’ 

%A/2  =  0, 

dx2  OK 

resulting  in  quadratic  distribution  of  temperature  and  potential 
T(x)  =  —  J- — x2  +  k\X  +  k2. 


(22) 

(23) 

(24) 

(25) 


and  a  constant  current  density 


J  = 


i 


with 

fcljpp^KBB,+n  k2=r(0))  fc3=,(0)> 

where  T(x)  and  tp(x)  atx  =  0 ,L  are  specified  as  the  boundary  conditions.  The  heat  flux  is  then  calculated  as 

Jq(X)  =  -fa/  +  +  «Jk^  X-Kk ,  +  ajk2.  (28) 

For  bilayered  thermoelectric  composites  consisting  of  two  different  materials  with  distinct  thermoelectric  properties, 
analytic  solutions  have  also  been  obtained  (Yang  et  al.,  2012).  In  this  work,  however,  we  are  interested  in  layered  composites 
with  periodic  microstructure,  which  we  approach  using  asymptotic  analysis,  as  we  discuss  next. 

3.  Asymptotic  homogenization 

We  consider  a  ID  composite  consisting  of  two  distinct  phases  layered  periodically,  as  schematically  shown  in  Fig.  1.  For 
such  a  composite,  two  different  length  scales  can  be  identified,  one  is  L,  the  macroscopic  length  of  the  composite  associated 
with  the  macroscopic  coordinate  x.  The  other  is  jy,  the  characteristic  length  of  the  composite  unit  cell,  for  which  a 
microscopic  coordinate  f  =  x/^  can  be  introduced.  While  the  material  properties  a(|),  a-(|),  and  k(£)  vary  fast  on  the 
microscopic  scale  periodically,  the  field  variables  such  as  T(x,f),  <P(x,£),  and  J0(x,  f)  vary  both  fast  on  microscopic  scale  and 
slowly  on  macroscopic  scale,  as  schematically  shown  in  Fig.  1(d).  What  we  are  interested  in  is  the  macroscopic  variation  of 
these  field  variables,  for  which  the  fast  microscopic  fluctuation  is  averaged  out,  so  that  the  effective  behavior  of  the 
composite  can  be  deduced. 

To  this  end,  we  expand  the  temperature,  potential,  and  heat  flux  into  polynomials  of  o 


T(X,  |)  =  T0(X,  I)  +  rjT ,  (X,  I)  +  n2T2(X,  f)  + 

<Kx,  I)  =  Mx,  f)  +  n4>  i(x,  0  +  v2Mx,  f)  + 

Jq(X,  I)  =Jqo(X,  I)  +  ^Qlix,  i)  +  rfjQ 2(x,  f)  +  ••• 

Inserting  these  into  governing  Eqs.  (19)  and  (21)  for  T  and  <-/>,  and  noticing  that 
A  ii* d 

dx  dx  ^  t/  <3|  ’ 

we  obtain  a  series  of  equations  grouped  by  the  orders  of  rj  as  following: 

dx  {K^fJ  +  dk  [~aJT°  +  K^+  KWJ  ~  W1 


(29) 

(30) 

(31) 


(33) 

(34) 


(35) 
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7  •  -J=a[^r+^)+aa\i^+^)’  (36) 

v„,v  1%JJ- dx  V ~'*J±0^''l)x^''l)l[)  V -“J'1  (37) 

We  first  examine  T0(x,£)  and  tf/0(x,  f),  the  lowest  order  of  field  variables.  By  integrating  Eq.  (33)  and  rearranging  the 
equation,  we  obtain 
dT0  _  C  i  (x) 

W~  W 

Since  T0(x,  |)  is  periodic  in  f,  and  the  integration  of  <JT0/<J|  over  the  unit  cell  vanishes,  it  is  concluded  that  Cj(x)  =  0,  and  thus 
T0(x,a  =  fo(x).  (39) 

Combining  this  with  Eq.  (34),  it  is  also  evident  that 

0o(x,|)  =  0o(x).  (40) 

In  other  words,  T0(x)  and  <p0(x)  do  not  vary  on  microscopic  scale,  and  they  represent  the  macroscopic  distributions  of 
temperature  and  potential  when  rj->  0.  They  thus  describe  the  effective  behavior  of  the  layered  composite  with  local 
fluctuations  of  temperature  and  potential  averaged  out,  which  we  seek  to  determine. 

We  then  examine  Ti(x,f)  and  (p\ (x,  S).  Making  use  of  Eqs.  (39)  and  (40),  we  can  integrate  Eq.  (35)  with  respect  to  |  as 


d£  ' 


"  dx  *(|)J 


x(l)  ’ 


(41) 


determined  again  from  the  periodicity  of  Tt(x,£)  on  the  unit  cell,  where  ( •  >  is  used  to  indicate  volume  averaged  quantities 
over  the  unit  cell.  Integrating  Eq.  (41)  one  more  time,  we  obtain 


Ti(x>s= Sr  hly1o-:h 

without  loss  of  generality,  si 


for  which  the  integration  constant  C3(x)  is  assumed  to  be 
Tq(x).  Similarly,  Eq.  (36)  can  be  recast  as 


which  by  integration  yields 


(45) 

As  such,  Eqs.  (43)  and  (45),  respectively,  give  T,(x,  |)  and  f/>,(x,  g)  that  fluctuate  fast  on  the  microscopic  scale  on  the  order  of  iy. 
This  will  be  sufficient  for  our  analysis,  and  higher  order  fluctuations  will  not  be  considered. 

In  order  to  derive  the  governing  equations  for  T0(x)  and  <j>0(x)  that  describe  the  effective  behavior  of  the  layered 
composite,  we  integrate  Eq.  (37)  over  the  unit  cell  to  obtain 

(46) 


which  can  be  simplified  by  using  Eqs.  (41),  (42),  and  (44),  resulting  in 

$+4K^2_(<f>2_<T>4>)-/2To=a  (47) 

This  is  the  governing  equation  for  the  macroscopic  temperature  distribution  T0(x)  in  the  layered  composite.  Furthermore,  we 
integrate  Eq.  (44)  over  the  unit  cell  to  obtain 

By  combining  Eqs.  (42)  and  (48),  we  can  rewrite  C2  as 

C2  =  -0~'(w+<Y}JTo  +  <a}J)’ 


(49) 
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which  can  be  used  to  recast  Eq.  (46)  a: 


This  is  the  governing  equation  for  the  macroscopic  potential  distribution  <j>0(x)  in  the  layered  composite. 

Eqs.  (47)  and  (50)  are  homogenized  governing  equations  for  temperature  and  potential  distribution  in  layered  periodic 
composite,  one  of  the  key  results  in  this  work,  and  the  coefficients  in  these  equations  can  be  easily  determined  from  the 
volume  averages  of  appropriate  quantities  in  the  unit  cell.  These  homogenized  governing  equations  are  considerably  more 
complicated  than  those  of  homogeneous  materials,  Eqs.  (22)  and  (23),  a  consequence  of  nonlinearly  coupled  transport 
equations  (Agoras  et  al.,  2009;  Ghosh  et  al.,  2009),  and  it  is  easy  to  verify  that  under  limiting  case  of  homogeneous  materials, 
the  homogeneous  governing  equations  are  recovered.  With  T0(x)  and  <p0(x)  solved  from  Eqs.  (47)  and  (50),  the  macroscopic 
distribution  of  heat  flux  can  be  derived  as 


iao  =  (^r1(<“W-dtJ),  (51) 

where  Eqs.  (41)  and  (42)  have  been  used.  Thus  the  problem  is  completely  solved,  if  we  can  determine  the  macroscopic 
temperature  and  potential  distributions  along  with  the  constant  current  density,  which  we  seek  to  solve  in  the  next  section. 
We  can  further  determine  the  local  fluctuations  of  temperature  and  potential  from  Eqs.  (43)  and  (45),  though  we  are  more 
interested  in  the  effective  behavior,  and  this  line  of  investigation  will  not  be  pursued  here. 


4.  The  effective  behavior 


4.1.  Field  analysis 

In  order  to  solve  for  the  temperature  distribution  in  the  layered  composite,  we  recast  Eq.  (47),  a  second  order  ordinary 
differential  equation  with  constant  coefficients,  as 

^+d1T0  +  d2  =  0,  (52) 

dx 

with 


d2={l){l}j2, 

which  can  be  solved  as 


Tq(x)  =  q,  cosiVcfx)  +  q2  sin(Vdix)- 


da 

di’ 


(53) 


where  q,  and  q2  are  constants  that  can  be  determined  by  boundary  conditions  as 
n  d'2  i  Tint  n  dz/ch  +  C0S(7d7L) 

<h  =j-+UU),  V2  =  - .  ,  rr-,. - ■ 

“i  sin(Vdil) 


(54) 


Not  surprisingly,  the  functional  variation  of  macroscopic  temperature  distribution  of  the  layered  composite  is  drastically 
different  from  the  quadratic  distribution  of  homogeneous  materials  in  Eq.  (24).  Having  derived  T0(x),  we  can  solve  for  <p0(x) 
from  Eq.  (48)  in  a  similar  manner  as 


(/>0  =  +  h2  sinfv^x)  +  h3  cos (\Zdjx)  +  h4, 


(55) 


with 


/A  = -<  ^  >;- (<  ^  r1  <  ^>2-<  ^  , 

h2  =  -( l  r1  < "  >g2  +  (<  l  r1  ( “  >2-<  ^ ^ 

h3=-(^>-1<“>qi-((^>-1(“>2-<^>)j^=, 


h4  =  <fl(0)-h3. 
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Note  that  both  temperature  and  potential  distributions  are  given  in  terms  of  yet  to  be  determined  current  density  J,  which 
can  be  solved  from  potential  boundary  condition  at  x=L, 

<j>0(L)  =  h\L  +  h2  sin (y/d^L)  +  h3  cos(\/di L)  +  h4.  (56) 

With  the  current  density  J  determined  as  such,  the  heat  and  energy  fluxes  can  be  determined  accordingly, 

Jqo  =  [ (-<Ji  Vdi~( " )J<?2) sin(yd7x)+(q2 757 -< ")/?,)  cos(y^x)  +  (57) 

Ju  =  <^r1(<f>rQ/-^)+^J.  (58) 

This  set  of  equations  completely  solved  the  macroscopic  variation  of  temperature,  potential,  current  density,  heat  flux,  and 
energy  flux  in  a  layered  composite,  and  thus  completely  describe  its  effective  behavior. 


4.2.  The  effective  properties 


The  effective  behavior  of  a  composite  material  can  often  be  described  in  terms  of  its  effective  constitutive  constants,  and 
this  can  also  be  attempted  for  thermoelectric  composite.  Since  the  effective  properties  of  a  nonlinear  composite  depends  on 
the  given  boundary  condition,  we  define  the  effective  thermoelectric  properties  through  the  following  equivalency  principle 
(Yang  et  al.,  2012).  Given  the  identical  boundary  conditions  of  temperature  and  electric  potential,  a  thermoelectric 
composite  with  a  set  of  effective  thermoelectric  properties  should  have  identical  current  density  and  energy  flux  as  a 
homogeneous  thermoelectric  with  the  same  set  of  properties.  With  such  equivalency,  it  is  clear  that  the  composite  and 
homogeneous  thermoelectrics  can  be  exchanged  under  the  specified  boundary  conditions. 

With  such  equivalency  principle,  we  examine  the  effective  electric  conductivity  first.  Consider  a  boundary  condition  of 
imposed  electric  potential  difference  only  with  AT  =  0,  and  compare  the  current  density  between  homogeneous  thermo¬ 
electric  and  periodic  composite,  we  conclude  that  the  effective  electric  conductivity  of  the  periodic  composite  is  given  by 


<7*(A0,AT=O)  = 


J 

A  <t>:0 


(59) 


with  the  current  density  determined  by  Eq.  (56).  Since  AT  =  0,  the  current  density  can  be  solved  analytically,  resulting  in 


(60) 


v=V<T><7>-<7>2- 


Clearly,  the  effective  electric  conductivity  depends  on  the  imposed  potential  difference  at  boundary  in  addition  to  the 
material  constants  of  the  constituents  and  the  volume  fraction/,  a  characteristic  distinct  from  linear  medium.  On  the  other 
hand,  if  we  impose  open-circuit  boundary  condition  where  J=0,  such  that 


(61) 


the  effective  Seebeck  coefficient  can  be  derived  as 


In  addition,  under  the  condition  of  J=0,  we  have 

,  _  /S-idTo 


from  which  the  effective  thermal  conductivity  can  be  derived  as 

r* _  Jqo  _  / 1  \-i 

[T(L)-T(0)]/L  V  ' 


(62) 


(63) 


(64) 


Notice  that  this  effective  thermal  conductivity  is  identical  to  that  of  a  layered  linear  thermal  conductor  uncoupled  from 
electric  conduction,  since  open  circuit  condition  is  assumed,  and  different  effective  Seebeck  coefficient  and  thermal 
conductivity  will  be  resulted  if  different  electric  boundary  condition  is  imposed.  From  these  effective  thermoelectric 
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properties,  the  effective  thermoelectric  figure  of  merit  is  then  defined  as 


Z*  = 


o*a*2 

K* 


(65) 


though  its  connection  with  the  thermoelectric  conversion  efficiency  remains  to  be  examined,  which  we  analyze  in  the  next 
subsection. 


4.3.  Thermoelectric  conversion  efficiency 


While  the  figure  of  merit  is  directly  related  to  thermoelectric  conversion  efficiency  for  a  homogeneous  material,  our 
analysis  on  bilayered  thermoelectric  composite  indicates  that  it  is  not  well  defined  for  heterogeneous  thermoelectric.  As  a 
result,  we  analyze  the  thermoelectric  conversion  efficiency  of  the  periodic  layered  composite  directly,  based  on  the 
thermoelectric  module  schematically  shown  in  Fig.  2.  One  leg  of  the  module  is  made  of  material  of  interest,  which  is  the 
layered  composite  with  periodic  microstructure,  and  the  other  leg  is  made  of  a  fictitious  material  with  zero  Seebeck 
coefficient  and  thermal  conductivity,  yet  infinite  electric  conductivity,  and  thus  it  only  serves  as  a  path  for  electric  current, 
not  involved  in  energy  conversion.  The  conversion  efficiency  of  this  idealized  thermoelectric  module,  as  a  result,  measures 
only  the  performance  of  the  composite  material  of  interest,  and  it  can  be  defined  as 


h_(JA)2R_  (JA)2R 
Eo  A/ulx  —  o  ’ 


(66) 


where  the  numerator  denotes  the  electric  energy  delivered  to  the  load  resistance,  whereas  the  denominator  denotes  the 
heat  (energy)  flowing  in  at  the  joint  from  the  hot  reservoir  with  temperature  T( 0).  From  Eq.  (58),  it  is  derived  that 

Ju\x  =  o  =  ( \  J"1  [( “  )T(0)J~  ^  I*  =  o] ,  (67) 

and  the  current  density  J  can  be  determined  from  integral  of  V<f>  around  the  circuit, 


(68) 


which  replaces  the  imposed  potential  at  boundaries  we  used  earlier.  The  efficiency  can  be  optimized  with  respect  to  the  load 
resistance  R,  and  for  homogeneous  materials,  this  lead  to  classical  equation  relating  figure  of  merit  ZT  to  optimal  conversion 
efficiency  under  given  temperature  difference, 


(69) 


where  Th  and  Tc  refer  to  temperatures  at  hot  and  cold  ends,  respectively.  For  composite  materials,  however,  such  simple 
relationship  is  no  longer  available,  and  we  have  to  evaluate  the  optimal  conversion  efficiency  numerically  in  general. 
Alternatively,  it  is  sometimes  insightful  to  recast  the  conversion  efficiency  as 


H  mj 


(70) 
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m =-<l  )Jl-(  “><  Ir1  [T(L)-m]  +  (V  r’<“  >2-<  ~  >);  Jql  t0  dx,  (7i) 

where  </;(0)  =  0  is  assumed  without  loss  of  generality.  The  efficiency  can  then  be  optimized  equivalently  with  respect  to  the 
current  density  instead,  though  the  current  density  is  upper  bounded  by  the  current  density  achieved  without  the  load 
resistance.  These  two  approaches  for  deriving  the  optimal  conversion  efficiency  are  equivalent,  and  both  requires  numerical 
computation  in  general,  yet  the  one  using  current  density  is  more  efficient  numerically. 

5.  Results  and  discussions 

5.1.  Field  distributions 

To  demonstrate  the  analysis,  we  consider  a  periodic  thermoelectric  composite  consisting  of  Bi2Te3  and 
Ag(Pbi_ySny)mSbTe2+m,  with  their  thermoelectric  properties  listed  in  Table  1  and  the  volume  fraction  of  Bi2Te3  being  /. 


Thermoelectric  properties  of  Bi2Te3  (Antonova  and  Looman,  2005)  and  Ag(Pbi_ySny)mSbTe2+m  (Androulakis  et  at,  2006). 


Material 

a  (xlO-6  V/K) 

a  (xlO3  S/m) 

*  (W/m/K) 

Bi2Te3 

200 

no 

1.6 

Ag(Pb,_j,Sny)mSbTe2+m 

270 

22 

0.77 

Fig.  3.  The  distributions  of  (a)  temperature,  (b)  electric  potential,  and  (c)  heat  flux  in  thermoelectric  composite  under  an  imposed  temperature  difference 
of  T(0)  =  300  K  and  T(L)  =  1 000  K,  with  ^(0)  =  0(1)  =  0  and  /  =  0, 0.4, 1 . 
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Notice  that  both  phases  have  excellent  Seebeck  coefficient  that  are  comparable  to  each  other,  yet  Bi2Te3  has  relatively 
high  thermal  conductivity,  while  Ag(Pbl_ySny)mSbTe2+m  has  relatively  low  electric  conductivity,  both  of  which  are  not 
desirable  for  high  thermoelectric  conversion  efficiency.  Consider  either  only  a  temperature  difference  or  an  electric 
potential  difference  is  imposed,  and  the  corresponding  distributions  of  temperature,  electric  potential,  and  heat  flux 
are  shown  in  Figs.  3  and  4,  where  /=  0.4  is  assumed,  and  homogeneous  materials  with/=0, 1  are  also  included  for 
comparison.  It  is  observed  that  under  both  types  of  boundary  conditions,  the  field  distributions  are  highly  nonlinear,  as 
expected,  and  the  heat  flux  is  not  constant.  Notable  temperature  rising  is  also  observed  in  Fig.  4  inside  of  the 
composite  under  imposed  potential  difference,  due  to  Joule  heating,  and  the  increase  in  temperature  is  most 
prominent  in  Bi2Te3,  which  has  higher  electric  conductivity,  and  its  heat  flux  is  also  larger  due  to  higher  thermal 
conductivity. 

It  is  worth  noting  that  the  quantitative  variations  of  temperature  and  electric  potential  for  the  layered  composite  are  not 
much  different  from  those  of  homogeneous  materials,  despite  drastically  different  functional  forms.  This  motivates  us  to  fit 
the  distributions  of  temperature,  potential,  or  the  heat  flux  in  the  composite  using  homogeneous  solutions,  Eqs.  (24),  (25), 
and  (28).  The  resulting  sets  of  thermoelectric  properties  identified  from  these  fittings  are  listed  in  Table  2,  where  it  is 
observed  that  different  curves  resulting  in  thermoelectric  constants  that  differ  in  orders  of  magnitude,  for  both  types  of 
boundary  conditions.  In  other  words,  it  is  difficult  to  identify  a  consistent  set  of  effective  thermoelectric  properties  that  fit 
the  distributions  of  temperature,  potential,  and  heat  flux  simultaneously  under  given  boundary  conditions.  The  comparison 
of  original  distributions  and  fitted  results  shown  in  Figs.  5  and  6  also  confirm  this  observation.  For  example,  in  Fig.  5,  it  is 
observed  that  the  data  fitted  by  heat  flux  leads  to  good  agreement  for  temperature  distribution,  yet  results  in  poor 
agreement  for  potential.  In  Fig.  6,  on  the  other  hand,  the  data  fitted  by  heat  flux  results  in  poor  agreement  for  both 
temperature  and  potential.  Similar  conclusion  can  be  drawn  for  data  fitted  by  temperature  or  potential.  This  highlights  the 
difficulty  in  homogenizing  the  effective  thermoelectric  properties  of  the  composite  that  would  resemble  the  behavior  of 
homogeneous  materials. 
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5.2.  The  effective  properties  and  conversion  efficiency 

Since  it  is  impossible  to  find  a  set  of  effective  thermoelectric  properties  of  composite  that  will  fit  its  distributions  of 
temperature,  potential,  and  heat  flux  into  homogeneous  solution,  we  calculate  the  effective  properties  of  the  composite 
versus  the  volume  fraction  of  Bi2Te3  using  equivalency  principle  instead,  as  shown  in  Fig.  7,  which  ensures  the  equivalency 
in  current  density  and  energy  flux  between  composite  and  homogeneous  material  with  the  identical  set  of  properties.  The 
effective  Seebeck  coefficient  and  thermal  conductivity  are  calculated  under  open  circuit  condition,  and  the  effective  electric 
conductivity  is  evaluated  with  T(0)  =  T(L)  =  300  K  and  Arp  =  0.2, 0.35, 0.5, 1 .0  K.  It  is  observed  that  there  is  a  noticeable 


Table  2 

By  fitting  curves  of  temperature,  electric  potential,  and  heat  flux  of  composite 
to  homogeneous  results,  sets  of  thermoelectric  properties  are  derived. 


Curves  fitted  a  ( x  1 0-6  V/K) 

a  (xlO4  S/m) 

X  (W/m/K) 

Under  imposed  AT : 

T0  541 

635 

890 

Po  467 

440 

838 

lost  29  077 

0.0002 

1 

Under  imposed  A 0: 

T0  2424 

3156 

957 

0o  1016 

809 

985 

!qo  276 

3 

620 

— Original 
T,  fitted 

■  «0  fitted 

JM  fitted 

0  2  4  6  8  10 


x  (mm) 


Fig.  5.  Distribution  of  (a)  temperature,  (b)  electric  potential,  and  (c)  heat  flux  of  thermoelectric  composite  are  fitted  to  the 
solutions  of  homogeneous  material  for  T(0)  =  300  K,  T(L)  =  1000  K,  0(0)  =  <p(L)  =  0,  and  /=  0.4. 
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x  (mm) 

Fig.  6.  Distribution  of  (a)  temperature,  (b)  electric  potential,  and  (c)  heat  flux  of  thermoelectric  composite  are  fitted  to  the 
solutions  of  homogeneous  material  for  /-/?(())  —  0,  ^(I)  =  0.1  V,  T(0)  =  T(L)  =  300  K,  and/=0.4. 


dependence  of  the  effective  electric  conductivity  and  figure  of  merit  on  boundary  condition.  We  further  calculated  the 
optimized  conversion  efficiency  of  the  composite  versus  the  volume  fraction  of  Bi2Te3  with  T(0)  =  800K  and 
T(t)  =  300  K,  as  shown  in  Fig.  8,  along  with  the  conversion  efficiency  calculated  using  the  effective  figure  of  merit 
based  on  Eq.  (69),  and  it  is  clear  that  there  is  a  large  discrepancy  between  the  actual  efficiency  and  the  efficiency 
evaluated  from  the  effective  figure  of  merit,  which  depends  on  boundary  condition.  This  suggests  that  the  effective 
thermoelectric  figure  of  merit  of  composite  is  not  well  defined,  and  does  not  correlate  with  its  thermoelectric 
conversion  efficiency. 

The  lack  of  correlation  between  the  effective  figure  of  merit  and  the  thermoelectric  conversion  efficiency  can  be 
understood  as  following.  To  define  the  effective  Seebeck  coefficient,  electric  conductivity,  and  thermal  conductivity,  two  sets 
of  boundary  conditions  are  required,  making  it  impossible  to  match  the  actual  optimal  boundary  condition  under  which  the 
conversion  efficiency  is  evaluated.  This  motivates  us  to  suggest  an  alternative  way  to  define  the  effective  thermoelectric 
properties,  using  equivalency  principle  under  optimal  working  condition  in  combination  with  the  requirements  of  Eqs.  (65) 
and  (69)  instead,  and  interestingly,  this  leads  to  two  sets  of  effective  properties,  as  shown  in  Fig.  9  along  with  the  effective 
Z*.  This  again  suggests  that  the  effective  properties  of  thermoelectric  composites  are  ill-defined. 


6.  Concluding  remarks 

In  conclusion,  we  developed  rigorous  nonlinear  asymptotic  homogenization  theory  to  analyze  the  coupled  transport  of 
electricity  and  heat  in  thermoelectric  composite  materials,  with  which  we  solved  for  the  macroscopic  field  distributions  that 
are  drastically  different  from  those  in  homogeneous  materials,  analyzed  the  overall  conversion  efficiency  using  an  idealized 
thermoelectric  module,  and  showed  that  the  effective  thermoelectric  properties  are  ill-defined,  and  the  effective  thermo¬ 
electric  figure  of  merit  does  not  directly  correlate  with  the  thermoelectric  conversion  efficiency.  The  analysis  thus  sheds 
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Fig.  8.  Conversion  efficiencies  of  composite  versus  volume  fraction.  There  is  a  discrepancy  between  directly  optimized  value  and  the  results  converted 
from  effective  figures  of  merit  from  different  boundary  conditions. 


considerable  insight  into  the  effective  behavior  of  thermoelectric  composites  that  can  be  used  for  their  design  and 
optimization. 
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